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Abstract. We investigate gravitational microlensing of 
point-like lenses surrounded by diffuse gas clouds. Besides 
gravitational bending, one must also consider refraction 
and absorption phenomena. According to the cloud den- 
sity, the light curves may suffer small to large deviations 
from Paczyhski curves, up to complete eclipses. Moreover, 
the presence of the cloud endows this type of microlensing 
events with a high chromaticity and absorption lines rec- 
ognizable by spectral analysis. It is possible that these ob- 
jects populate the halo of our galaxy, giving a conspicuous 
contribution to the fraction of the baryonic dark matter. 
The required features for the extension and the mass of 
the cloud to provide appreciable signatures are also met 
by several astrophysical objects. 

Key words: Gravitational Lensing, Gas Clouds, Circum- 
Stellar Matter, Pre-Main Sequence Stars, Planetary Neb- 
ulae, AGB and post-AGB Stars 



1. Introduction 

In the last years gravitational microlensing has rapidly 
grown to an important astrophysical tool with wide ap- 
plications. In particular, it gave some evidence for the 
presence of baryonic dark matter, in form of compact ob- 
jects (MACHOs), in the Galatic halo (Alcock et al. 2000a, 
Lasserre et al. 2000) although the number of microlensing 
events in the direction of the Magellanic Clouds, reported 
by the observational teams MACHO and EROS, is lower 
by a factor ~ 5 as predicted by the standard halo mod- 
els. The estimated MACHO fraction is now about 20% 
of the total Galactic dark matter mass, though it must 
be kept in mind that this value is still affected by a very 
high uncertainty. On the other hand several hundred of 
events have been found towards the Galactic bulge region 
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(Alcock et al. 2000b, Wozniak et al. 2001) and few events 
towards some spiral arm regions (Derue et al. 2001). Most 
likely the lenses are in these cases low mass stars or brown 
dwarfs. 

Classical microlensing consists of the magnification of 
the luminous flux coming from a star due to the pres- 
ence of a deflecting mass close to the line of sight. Within 
this picture, microlensing yields an indirect observational 
method for the detection of non-luminous bodies. 

However, compact objects are not the only candidates 
for baryonic dark matter in the halo. Some years ago, 
it has been proposed that an appreciable fraction of the 
Galactic dark matter might be in form of self-gravitating 
gas clouds with masses of the order of Jupiter's mass 
(10~ 3 M Q ) and radii R w 10 AU (Pfenninger et al. 1994; 
De Paolis et al. 1995a, 1995b, 1996, 1998a, 1998b; Gerhard 
& Silk 1996). For diffuse objects two more phenomena 
must be taken into account. Traveling through a gaseous 
medium, light basically undergoes refraction and absorp- 
tion. So, even if these gas clouds would produce no appre- 
ciable microlensing effect, classical refraction could sup- 
ply an alternative amplification mechanism. The resulting 
light curves resemble gravitational microlensing ones, but 
can be distinguished by their peculiar features (Drainc 
1998; Rafikov & Draine 2001). 

In this paper, we mean to study an intermediate pos- 
sibility, which might be of relevance for the estimate of 
baryonic dark matter in the halo. After almost ten years 
of microlensing campaigns, at a good confidence level, we 
can say that MACHOs exist. If gas clouds are effectively 
present in the halo, they could be likely associated to MA- 
CHOs. If this is the case, these Pointlike Lenses Associated 
to Gas (hereafter, PLAG) would be missed by standard 
microlensing observations, because their light curve would 
be significantly altered by the surrounding gas cloud. On 
the other hand, isolated gas clouds could be quite difficult 
to detect without a powerful amplification mechanism. 

In PLAGs gravitational lensing by the central object 
would provide this mechanism, but the secondary effects 
due to the cloud would show up in a weak or strong fash- 
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ion, depending on the density of the diffuse matter. Clas- 
sical absorption by Rayleigh scattering would essentially 
block all light attempting to cross high density clouds. 
For lower densities, a moderate absorption would be con- 
trasted by microlensing and refraction both concurring in 
the amplification of the flux coming from the source. In 
this case, a considerable reddening of the light must be ex- 
pected, producing a high chromaticity in the light curve. 
Of course, the existence of the gas cloud could be revealed 
without doubts by taking a spectrum of the source during 
the microlensing event. In this way, the chemical composi- 
tion of the matter surrounding the central compact object 
can be determined. 

Several studies about the stability of gas clouds seem 
to favour the formation of a steady state in protostellar 
nebulae (Cazes & Tohline 1998). If MACHOs were formed 
by the collapse of gas clouds, a large portion of this cold 
gas may still surround some of them, with a reasonable 
formation of PLAGs. 

The abundance of dark PLAGs in the halo is not 
known, leaving open several possibilities. They could be 
completely absent, if some instability mechanism inter- 
venes. They may be a negligible fraction of MACHOs, if 
their mean life is too small. On the contrary, they could be 
relatively abundant and partly explain the discrepancy be- 
tween observations and theoretical predictions about the 
microlensing rate, since the missing microlensing events 
would be actually PL AG events, not detected by usual 
observations. If this is true, then microlensing experiences 
should be properly modified in order to detect PLAG 
events efficiently. 

Besides dark PLAGs in the halo, many well-known as- 
trophysical compact objects surrounded by gas are present 
in the galactic disk. We can mention, for instance, plane- 
tary and proto-planetary nebulae, proto-stars, LBV neb- 
ulae, nebular variables (T Tau stars and R W Aur stars), 
eruptive variables (R CrB stars, Wolf-Rayet stars, novae), 
and so on. Of course, not all of them are able to produce 
microlensing effects that are really observable by monitor- 
ing programs, since their parameters should satisfy several 
constraints. 

A fundamental tuning is required between the size of 
the cloud R c and the Einstein radius Re- The ratio R c /Re 
should be in the range I _1 -j-10 2 . This tuning is demanded 
by the necessity of having a microlensing time comparable 
to the time required by the source to pass behind the 
whole gas cloud. In fact, if the cloud is too small, it would 
act, at all effects, as a compact object. On the contrary, 
if the cloud is much larger than the Einstein radius, with 
the typical velocities of galactic objects, the time taken 
for the source to pass behind the whole cloud may easily 
rise up to years. However, even in this case, these long 
duration events could be included among the objectives 
of observational campaigns following stars for very long 
periods. 



For R c /Re > I0 2 , we would observe no deviation from 
the usual microlensing light curve, since the effects of the 
cloud would be manifest after decades and could not be 
practically noticed. No information on the size and the 
shape of the cloud would be retrieved in this case. 

In Sect. 2, we introduce the basics of the three phenom- 
ena affecting the light passing through the cloud: gravita- 
tional lensing, refraction and absorption by Rayleigh scat- 
tering. In Sect. 3, we analyze the general features of the 
light curves expected from PLAGs. In Sect. 4, we briefly 
discuss the properties of the afore-mentioned PLAG can- 
didates. Sect. 5 contains the summary. 

2. Light passing through gas clouds: basic 
phenomena 

We restrict our analysis to spherical gas clouds surround- 
ing point-like masses. This simplification allows an exten- 
sive analytical investigation. Moreover, spherically sym- 
metric PLAGs should be easier to distinguish from vari- 
able stars and other usual background to microlensing. 

The lens equation contains both refraction and gravi- 
tational lensing. The latter can be separated into two con- 
tributions: the deviation due to the central object and the 
one produced by the gas cloud. Thanks to the spherical 
symmetry, we can write the one dimensional equation 

y = x- D °i Dls [a (x) + a c {x) + a r (x)] , (I) 

J-Jos 

where 

, \ 4GM 

«<>(*) = — (2) 

is the deflection angle produced by the gravitational field 
of the compact object (of mass M) at the center of the 
PLAG (Schneider, Ehlers and Falco 1992); a c (x) is the de- 
viation angle due to the gravitational field of the gas cloud; 
a r {x) is the contribution to the deviation coming from re- 
fraction of light passing through the cloud. As usual, y 
(x, respectively) is the distance of the source (of the im- 
age, respectively) to the optical axis, defined as the line 
connecting the observer and the lens. 

The splitting of the deviation angle into the sum of 
three separate contributions is possible in the linear ap- 
proximation. Higher orders would mix up the three effects. 

The gas cloud is described by a continuous density 
distribution p(r) which is a function of the radial distance 
from the center of the PLAG, according to the spherical 
symmetry hypothesis. Since the PLAG has limited size, 
we call R c the value of the radius where the density falls 
to zero. The total mass of the cloud is M c . The projected 
density is 

£(&)=/ p(Vx 2 + z 2 )dz. (3) 
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The gravitational bending for light rays passing out- 
side the cloud (x > R c ) is just 4 ^f c , while for x < R c , a c 
decreases to zero. It is generally given by 



ac (x) = A 4^[\(x')M = ^ 
cr x J c z x 



(4) 



where m(x) is the mass of the portion of the cloud con- 
tained in the cylinder of radius x (Schneider, Ehlers and 
Falco 1992). 

Light rays passing inside the gas cloud undergo 
Rayleigh scattering by neutral gas (essentially Hydro- 
genum and Helium) and eventually dust. The macroscopic 
consequences are refraction and extinction, both depend- 
ing on the density of the neutral gas /?„(r) and on its 
chemical composition. 

The refraction index of a gas is 



n\(r) = 1 + a\p n (r), 



(5) 



where a\ is the specific refractivity. The deviation angle 
can be easily derived by the Fermat principle assuming 
the deviation to be small. Its expression is 



a r (x) 



-ax- 



dTi„ 
dx 



= —2axx 



dr dp n 



2 dr 



r'- — x 



(6) 



where £„ is the projected density of neutral gas (confront 
Draine 1998, Eq.l). 

The magnification of the images follows from the one 
dimensional lens equation ([!]) 



(7) 



with y given by the RHS of Eq. ([!]). This magnification 
factor is a consequence of the geometrical bending on the 
light rays, caused by refraction and gravitational field. 

The observed intensity is moreover modified by com- 
bined effects of scattering and absorption, giving rise to 
the so-called extinction, described by the Lambert's law 



h = Ix,oe~ 



(8) 



where /a,o is the intensity that we see in absence of ab- 
sorption along the line of sight, i.e. before and after the mi- 
crolensing event, and Tx is the optical depth of the cloud, 
given by 



(9) 



where kx is the absorption coefficient. 

Summing up, each image will have a total amplification 
given by the geometrical magnification (R) due to light 
bending (both gravitational and refractive) weighted by 
the transmission coefficient e~ TA . 



3. General properties of PLAGs 

The fundamental parameters of the PLAG are the mass 
of the central object M, the radius of the cloud R c and its 
mass M c . Most statements about the aspect of the light 
curves can be given in terms of these parameters, or their 
derived quantities, without need to specify the shape of 
the density profile (if it is smooth enough). 
We define the Einstein radius 



Re 



UG(M + M c ) DqlDls 



Dos 



(10) 



referring to the total mass of the system. 

In this section, using a simplified model, we will discuss 
the phenomenology that should be generally expected in 
microlensing events caused by PLAGs. 

3.1. Sample family of PLAGs 

If we consider PLAGs formed by hot nebulae, in the cen- 
tral region the gas maybe ionized, being uncffcctivc for 
classical Rayleigh scattering, suggesting p n < p. We are 
mainly interested in cold dark PLAGs; so, in a first ap- 
proach, we will assume that this region, if present, is small 
enough with respect to the total cloud and the Einstein 
radius, so that p n ~ p. In this way, we can rapidly obtain 
a simple description of the basic phenomenology. Never- 
theless, the extension to the case p n < p is straightforward 
and can be easily faced successively. 

The dependence on the nature of the gas is not so 
strong to substantially affect the discussion, since the or- 
der of magnitude of the specific refractivity and the ab- 
sorption coefficient is the same for all the interesting el- 
ements or molecules. Thus, in all practical examples, we 
will assume a mixture of H2-He with 24% He by mass. At 
A = 4400 A for this mixture, we have (AIP Handbook 
1972, following Draine 1998) 

1.243cm 3 g' 1 



ax 



5.4 x 10 



-5 



0.44/iTO\ o i 

r cm q 

A J 



(11) 
(12) 



The variety of density profiles in nebulae is very large 
and most of them have no satisfying description at the 
present time. This variety seems to make the choice of the 
density model a hard task. However, once projected along 
the line of sight, all realistic density distributions are gen- 
erally smoothly peaked at the center and decay towards 
the border in a more or less steep fashion. For this rea- 
son, the basic phenomenology of PLAG events does not 
depend on the specific shape of the original density pro- 
file but can be discussed by analyzing one representative 
model. 

We choose the function 







1 - 



for r < R c 
otherwise 



(13) 
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parameterized by the exponent f3 > 0, to describe the den- 
sity profile of a typical PLAG. The choice of this model 
is strongly motivated by the fact that most physically in- 
teresting quantities can be calculated analytically starting 
from this distribution, while all the general features of a 
PLAG microlensing event can be obtained. Moreover, the 
parameter (3 can be adjusted to suit to distinct physical 
situations: high /3's would produce clouds concentrated at 
the center with smooth boundaries which would appro- 
priately describe stable or nearly stable systems; low /3's, 
on the contrary, can be used for expanding clouds, which 
are generally characterized by a steep boundary and where 
abrupt changes in the density may produce discontinuities 
in the refraction angle (see Sect. 3.4 and Fig. |l|b). 

Of course, if one is interested in the details of the mi- 
crolensing light curves, it is necessary to know the real 
shape of the nebula; but for all the highlights of microlens- 
ing light curve, the study of this sample model will be 
sufficient. 

The central density can be expressed in terms of the 
mass of the cloud 



4nr pp{r)dr 



1] 



r[j9- 



The projected density is 

^/ 2 i? c r[/3 + i] 



r[j9- 



(14) 



(15) 



The gravitational bending due to this density distri- 
bution is 



a c (x) 



4Gp° g ir 3 / 2 R 3 c T [(3 + 1] 



2 xr [p - 



The refraction angle is 



a r {x) = ot\ 



2[3p 8 ^/ 2 R 3 c T\ 



T[J3- 



R, 



x 2 \ 



(16) 



(17) 



Of course, we cannot take as physical the divergence 
for < [3 < 1/2. The simple expression we are using 
for the refraction angle is valid for small deviations, so 
the presence of such divergence signals a breakdown of 
this hypothesis and the need for higher order corrections. 
These corrections should provide a cutoff for the refraction 
angle. 

In Fig. [I] we plot the gravitational bending, the refrac- 
tion angle and the transmission coefficient of our cloud 
model for M c = 0.01 M , R c = 6 AU and for different 
values of (3. 

3.2. Light curves 

We can distinguish two main classes of PL AGs, according 
to the ratio between R c and Re- 
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Fig. 1. (a) Gravitational bending due to the cloud, (b) Re- 
fraction angle, (c) Transmission coefficient. All quantities 
are plotted for different values of (3. 



In small PL AGs (R c < Re), the positive parity image 
(the image formed outside of the Einstein ring) is not af- 
fected by the existence of the gas cloud. Its position is just 
that of the principal image of a point-like lens with mass 
M + M c . Since light passes outside of the cloud, neither 
refraction nor absorption modify the simple gravitational 
lensing effect. 

The secondary image has negative parity and is formed 
inside the Einstein ring. If the source is close enough to 
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Fig. 2. The six categories of light curves expected for PLAG events as discussed in Sec. 3. 2. The dotted lines represent 
the classical Paczyhski light curves for gravitational lensing by the total mass M + M c . The common parameters, 
typical for a disk/bulge microlensing events, are [3 — 5/2, D i = 8kpc, Di s — 2kpc, M — 1 M , impact parameter 
b = {).2Re and transverse velocity vt = 200km/s. The other parameters are: (a) M c — O.OIMq, R c — 3.3 AU; (b) 
M c = 0.01M©, R c = 4.5 AU; (c) M c = O.O3M , R c = 6 AU; (d) M c = 0.7M Q , R c = 8 AU; (e) M c = 0.5M©, R c = 12 
AU; (f) M c = 0.2M Q , R c = 6 AU. 



the optical axis, this image is external to the cloud as well. 
In this situation, the PLAG acts as a point-like lens at all 
effects. When the source is far from the optical axis, the 
secondary image is formed inside the PLAG. In this case, 
its position is modified by the non trivial form of the de- 
viation angles a r and a c . Moreover, the light flux coming 
from this image is partially absorbed. However, the sec- 
ondary image is relevant for the total flux only when it is 
close to the Einstein ring, otherwise it is highly demag- 
nificd. So, these additional effects intervene in a regime 
where the image does not contribute much to the light 



curve. On the basis of these considerations, the typical 
light curve we expect for R c < Re is a Paczyhski- like 
curve slightly lowered far from the maximum, where the 
secondary image enters the cloud. This deviation can be 
easily confused in the background photometric noise. 

To ensure the effect to be appreciable, the size of the 
cloud should be as close as possible to the Einstein ring 
and the density as large as possible. In fact, in this case, 
the secondary image would remain inside the cloud during 
a large fraction of the microlensing event. A large absorp- 
tion would depress it until it exits from the cloud. So, 
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just the region in the neighborhood of the peak would be 
pure microlensing, while outside of this interval, the sec- 
ondary image would be depressed, changing the shape of 
the curve. An example is given in Fig. |2|a, drawn (like the 
others) using our sample model. 

The time of the entrance of the secondary image into 
the cloud is 



t = t + 



1 



Vt 



'3l 



b 2 



(18) 



where b is the impact parameter of the source to the opti- 
cal axis, Vt is the relative projected velocity and to is the 
time of the closest approach. 

An interesting phenomenology arises when R c > Re- 
The secondary image is always inside the cloud, being dis- 
placed and absorbed. The principal image would enter the 
cloud only in a neighborhood of to, for sufficiently small 
impact parameters. The time of the exit of the principal 
image from the cloud is still given by Eq. (|lj) . Depending 
on the density of the cloud, we can distinguish some cases. 

If the density is low, the absorption slightly lowers the 
maximum with respect to Paczyhski (Fig. |^b) . 

For moderate densities, the curve will be depressed 
when the principal image enters the cloud. So, during the 
ascent to the maximum, the curve will suddenly decrease 
its derivative and then complete its rise to a considerably 
smaller maximum (Fig. |^c). 

For high densities, the absorption will prevail on the 
magnification. When the principal image enters the cloud, 
the light reaching the observer falls rapidly. Eventually, 
the magnification at t = to can still produce a low cen- 
tral peak in the valley (Fig. ||d), but for higher densities, 
even this faint feature would disappear. If the time of the 
entrance is far from to t the event resolves in a complete 
eclipse of the source (Fig. ^e); otherwise, the initial part 
of the ascent can signal the presence of the central lens 
inside the cloud (Fig. ||f). 

The curves just described represent what can be gen- 
erally expected from this physical situation. As already 
quoted, the specific shape obviously depends on the den- 
sity profile considered. This fact can be exploited to try 
a reconstruction of the density profile from the details of 
the microlensing light curve. 

We can make some simple considerations about the 
relative weight of gravitational lensing, refraction and ab- 
sorption in the light curves, A(t). Keeping (3 — 5/2, we 
have compared the partial light curves (obtained neglect- 
ing one of the above-quoted phenomena at a time) with 
the exact ones. First, for each value of M c /M and R c /Re, 
we have calculated the fractional deviation of the curve 
without cloud gravitational lensing with respect to the 
complete one, in the course of the whole microlensing 
event 



/ac=o(*) 



\A(t) - A(t)\ ao=0 \ 
A(t) 



(19) 



In Tab. [I] we report the maximal deviation defined as 
max {fa =o(t)} occurring in each curve. In the same way, 
in Tab. ||, we have collected the maximal fractional devi- 
ation of the curves without cloud refraction a r with re- 
spect to the complete ones. Tab. [| contains the maximal 
fractional absorption. All light curves were computed as- 
suming an impact parameter b = 0.2Re, M = IMq and 
with A = 4400 A. 

With these tables it is easy to distinguish the regimes 
where the different phenomena become important, keeping 
the gravitational lensing of the central object as reference. 

Cloud gravitational lensing can be neglected when the 
cloud mass is below a few hundredth of the central mass. 
Its effects become heavier for higher cloud masses, maxi- 
mally for radii of the order R c /Re ~ 1.5. 

Refraction is always of the same order as cloud gravi- 
tational lensing, dominating only for small radii. It is not 
a surprise that the sum of the fractional deviations due to 
refraction and gravitational lensing is sometimes greater 
than one. In fact, their action is by no way additive, since 
they simultaneously concur in the shift of the images and 
modify the magnification formula in a non-linear way. 

All ranges of the cloud parameters where the refrac- 
tion contribution is non-negligible are also affected by a 
large absorption. High masses completely block the light 
as we can read in the last row of Tab. | Cloud masses 
around one tenth of the central mass still show good pos- 
itive deviations coming from light bending, maintaining a 
moderate absorption. It is in this mass range that PLAGs 
have the highest probability of being detected. 

3.3. Chromaticity 

Both refraction and absorption have a moderate depen- 
dence on the wavelength of the incoming wave. This means 
that a strong chromaticity should be expected for PLAG 
events, mostly caused by a selective absorption. In a minor 
extent, blending and refraction may slightly modify the 
color. Given the A -4 dependence of the absorption coeffi- 
cient, Eq.(|rj), the light filtered through the PLAG would 
be strongly reddened, as it is evident from Fig.^J, where we 
reported the light curve of the example treated in Fig.^c 
at different wavelenghts: A = 4400 A (blue), Fig.||a, and 
A = 7000 A (red), Fig.|b. For this reason, the achromatic- 
ity test of classical microlensing surveys would automati- 
cally reject PLAG events and should be relaxed in order 
to find them. Moreover, the detection of a PLAG would 
be helped choosing filters operating in opposite bands of 
the spectrum so as to put in evidence the real nature of 
the chromaticity in an optimal way. 

3.4- Boundary discontinuities 

In expanding shells produced by old stars, the density pro- 
file is generally characterized by a steep boundary, where 
the density suddenly changes from zero to a value signifi- 
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Table 1. Maximal fractional contribution to the complete light curve coming from the gravitational bending due to 
the cloud. 
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Table 2. Maximal fractional contribution to the complete light curve coming from refraction. 
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Table 3. Maximal fractional absorption. 



cant for refraction and absorption effects. Of course, this is 
not the case for cold dark PLAGs on which we are mainly 
focused in this article, which are supposed to be stable or 
long-lived systems; but a wide range of alternative PLAG 
candidates (see Sect. 4) require an appropriate description 
of boundary effects. 

For this reason, in this subsection we check whether ~ 2v / 2/(-R c ) 
interesting observable effects could be expected when one 
of the images enters an expanding steep-front cloud. 

First of all, let us try a perturbative resolution of the 
lens equation for x ~ R c - We let 



with f(r) a regular non vanishing function at r — R c . This 
expansion can be performed in general, not only for the 
sample family studied so far, but also for a generic density 
profile. 

The projected density will go as 



tic 



(23) 



Expanding all deviation angles to first order in e, we 



get 



y = R c 

and 



Re x 
Tt~ Sy 

tic 



(20) 



do 



R% M 
~R~ C M + M C 



Rc 

e 



R% M c j 
R c M + M C \ R, 



x = R c - e. (21) 
In the neighborhood of R c , pp(r) falls to zero as 



a r ~ a x 2V2f(R c ) [0 + 



1\ £ P-i 



2/ p' 3 " 

tic 



(24) 
(25) 
(26) 



Pfiir) = f(r) 1 - 



Rr 



(22) 



which can be taken seriously only for f3 > ^ . For lower val- 
ues, a cut-off from higher order terms should be expected. 
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Fig. 3. Chromaticity of the PLAG microlensing event, (a) 
A = 4400A(blue). (b) A = 7000A(red). The dotted lines 
are again the classical Paczyhski light curves, (c) Ratio 
between the two curves. 



The perturbed lens equation reads 



Ri 



(27) 



and the magnification for an image close to the boundary 
is 



R E 

m 



2e 
R r 



A 

R r 



Fig. 4. Details of light curves at the entrance of an im- 
age into the cloud. The dashed curve is the classical 
Paczyhski. (a) Plots for (3 = 3,5/2,2,3/2, 1. (b) Plots for 
0= 1/2, 1/4. 



Rl ( 2s 
l + -^f 1 + — 



Rl 



(28) 



Refraction yields /3-dependent terms to the lens equa- 
tion and to the magnification. It is easy to analyze the 
boundary effects that may arise for low values of (3. In 
Fig. 4, the light curve is shown for different values of /3 at 
the time of the entrance of an image inside the cloud. 

According to the value of (3, we can have different 
cases. 

1. /3 > 5/2: The refraction terms introduce subdominant 
corrections. The magnification is a continuous func- 
tion. 

2. (3 = 5/2: The first derivative of the magnification has a 
discontinuity, being higher for r < R c . The light curve 
suddenly decreases its derivative when an image enters 
the cloud. 

3. 3/2 < (3 < 5/2: The refraction term becomes dominant 
in the magnification. The derivative of the magnifica- 
tion goes to infinity in R c and the light curve has a 
continuous depression when an image comes into the 
cloud. 



V. Bozza et al.: Microlensing by Compact Objects associated to Gas Clouds 



9 



■5. 



(3 = 3/2: The refraction terms in the magnification do 
not vanish on the border. The light curve has a finite 
discontinuity 

1/2 < j3 < 3/2: The refraction term in the second 



7. 



square brackets of Eq.(28) diverges and the magnifi- 
cation falls off to zero. Then it rises from zero to the 
value it would have without refraction. 
(3 = 1/2: The lens equation has no solution until Sy 
becomes greater than A. This fact provokes a hole in 
the light curve at the entrance of the image into the 
cloud, illustrated in Fig. 4b for our sample model with 
(3 = 1/2. 

< (3 < 1/2: Again we have a hole in the light 
curve. Increasing Sy further, two values of e solve the 
equation. One approaches the original linear relation, 
behaving as the continuation of the original image. 
The other approaches zero, i.e. an additional image 
is formed very close to the border of the cloud. 
So, when an image enters the cloud, after a short in- 
terval where it temporarily disappears, it comes out 
doubled, with the additional image staying close to 
the border. The formation of this pair happens on a 
radial critical curvef] and the magnification diverges at 
the formation of the two images. 



These boundary features produced by refraction are 
a very interesting tool to investigate the shape and the 
nature of the cloud; but they last few hours and would 
be hard to detect. It must be recalled that for such steep 
variations the finite size of the source intervenes to smooth 
peaks and holes. Looking at the their typical time scales, 
these features would not be more elusive than planetary 
microlensing, after all. For very sharp density profiles, 
higher order corrections in the refraction formula should 
be taken into account for a correct description at very low 
(3. 

3.5. Additional critical curves 

In the previous subsection, we have seen that the struc- 
ture of the Jacobian can be predicted once we know the 
expansion of the density distribution near the boundary 
(p2|). In this way, we have found that refraction would 
produce radial critical curves very close to the boundary 
in densities with very steep fronts^]. Are these the only 
possible additional critical curves? 



1 We recall that the critical curves are defined by the van- 
ishing of the Jacobian determinant of the lens application. In 
axisymmetric lenses a critical curve is tangential if y(x) = 0, 
while is radial if =0. 

OX 

2 Actually, for /3 < 1/2 and R c > Re, between the additional 
radial critical curve and R c there is also a tangential critical 
curve, but this feature is certainly a product of the small angle 
approximation in a r , since it comes for very high refraction 
angles. 



In isolated clouds, without the central object, critical 
curves and caustics can be produced when the density is 
sufficiently high, as discussed by Draine (1998). However, 
when a central mass is present and the cloud mass is non- 
dominant with respect to the central mass, it is far more 
difficult to create critical curves as we shall discuss in the 
following. 

To have a tangential critical curve, we need the RHS 
of the lens equation (Q) to vanish. Usually the Einstein 
ring is the only critical curve of this class. It is difficult 
to have additional curves inside the Einstein ring, because 
we would need a very highly negative refraction angle to 
compensate the other terms. This kind of refraction angle 
could only come out in shell densities with a very steep 
internal front. Additional tangential critical curves can be 
created outside of the Einstein ring with very high pos- 
itive refraction angles. Such angles are beyond the reach 
of a linear theory and higher order terms would cut off 
the refraction angle excluding the formation of tangential 
critical curves. 

The situation is similar for radial critical curves. Here 
we have to kill the derivative of the RHS of (Q), i.e. we 
have to satisfy the equation 



.. D iDi s . /r . . 
1+ — - [c*a£ (x) + 



D os 

AG 



M + m(x) 



- 2nZ(x) ] = 0. 



(29) 



Without the refraction term, this condition reduces 
to the usual gravitational lensing one (Subramanian & 
Cowling 1986). Refraction helps to reach this condition 
but it becomes effective at high densities. In the range 
of masses we are interested into, additional images would 
be formed only in optically thick clouds so that no sig- 
nature of their existence would be visible on the light 
curves. On the other hand, decreasing the central mass 
and the mass of the cloud correspondingly, the role of re- 
fraction becomes prominent while the cloud becomes more 
and more transparent. Finally, going to the regime studied 
by Draine (1998) (no central mass and M c ~ 1O~ 3 M ), 
caustic-events start to be definitely observable. However, 
in PLAGs, the presence of a relevant central mass pushes 
the threshold given by Eq. (^) far beyond a high density 
regime, where the event just looks like a complete occul- 
tation (Fig. |e). 

Again, only in steep fronts it is possible to have radial 
critical curves as we saw in Sect. 3.4 for (3 < 1/2. 



4. Physical Nature of PLAGs 

In this section we examine the possible PLAG candidates, 
on the basis of the phenomenology risen up to now. We 
have seen what the microlensing light curves should look 
like and we have established an optimal range for the pa- 
rameters. 
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As regards cold PLAGs in the halo, i.e. MACHOs as- 
sociated to gas clouds, they would be ideal candidates, 
since their size could be about of the same order of the 
Einstein radius. The existence of such objects in the halo 
would have very interesting consequences in our knowl- 
edge about the gas distribution in the halo. It could be 
argued that a large fraction of cold baryonic mass in the 
halo can be found as gas surrounding MACHOs, since the 
small masses of these objects would favour the stability of 
the clouds. If this is the case, the density of PLAGs may 
be comparable to ordinary MACHOs, and several PLAG 
events per year can be expected by observing towards the 
Magellanic Clouds. 

Besides cold dark PLAGs in the halo, which surely rep- 
resent the most interesting possibility, in this Section we 
will examine known types of nebulae to look for additional 
PLAG candidates. 

Circumstellar disks are commonly found at IR wave- 
lenghts around young stars and around low mass and in- 
termediate stars that are in their evolutionary stage from 
red giants to white dwarfs. 

Proto-stars could be good candidates. In fact, they are 
isolated interstellar clouds undergoing gravitational col- 
lapse. The clouds are actually flattened to disks with cir- 
cular or elliptical projected shape according to the tilt 
angle with respect to the line of sight. Proto-stars can 
be found in star-forming regions, usually located in the 
disk and in the spiral arms of the Milky Way. Of course, 
to observe light curves as those shown in Fig.^, we are 
mainly interested in Proto-stars with 1 < M/M c < 10 
and R c ~ 1 - 10 2 AU. 

In the last decade, a good number of Proto-stars have 
been detected in our Galaxy (see for instance Beckwith et 
al. 1990, Hillenbrand et al. 1992, McCaugherean, O'Dell 
1996, Schneider et al. 1999, Padgett et al. 1999, Tuthill 
et al. 2001) and, very recently, the HST found also a very 
compact star-forming region in LMC (Heydary-Malayeri 
et al. 2001). Twin proto-planetary disks, each of them with 
R c ps 20 AU, have also been observed (Giovanardi et al. 
2000), suggesting the possibility of binary PLAG events. 

Other good candidates could be the proto-planetary 
nebulae (PPNe). They represent a transitory stage, be- 
tween the asymptotic giant branch (AGB) and planetary 
nebulae phases, in the evolution of intermediate and low- 
mass stars. This stage is characterized by an expand- 
ing circumstellar envelope of gas surrounding a central 
star (CPPN). In our phenomenological description, we 
can consider the PPNe as static systems, because the 
expansion velocity of the circumstellar envelope is v e ~ 
25 fcmsec" 1 (Osterbrock 1989). The typical mass of the 
CPPN is M 0.6 Mq and the mass of the circumstel- 
lar envelope is of the same order. Since they are expand- 
ing systems, it is certainly possible to observe PPNe with 
R c ~ Re- PPNe were only discovered and studied in the 
last twelve years, mostly thanks to IRAS and HST (Ueta 



et al. 2000). For a good overview about these objects, see 
Kwok 2000. 

As well as in our galaxy, where we know roughly 10 3 
planetary nebulae (PNe), intracluster PNe (IPNe) have 
also been detected thanks to the surveys towards the Virgo 
cluster (Mendez et al. 1997), M87 and its surrounding halo 
(Ciardulo et al. 1998, Feldemeier et al. 1998). 

Besides Proto-stars and PPNe, we can consider other 
kinds of astrophysical objects that are associated with gas 
clouds. For instance, recently the VLT discovered a small 
cone-shaped nebula near the the old neutron star RX 
J1856. 5-3754. This nebula extends for roughly 60 AU (van 
Kerkwijk 2001), and the neutron star is located very close 
to the top of the cone. Similar bow wave nebula have been 
found around other fast-moving supernova remnant, like 
IC 443 (Olbert et al. 2001), radio pulsars, like PSR B1706- 
44 (PSR J1709-4428) and PSR B1643-43 (PSR J1646- 
4346) (E. Giacani et al. 2001), nova system, like BZ Cam 
(Prinja et al. 2000), and could be good PLAG candidates. 

All these objects are located mainly in the bulge, in 
the disk and in spiral arms of galaxies. So, good targets 
for the search of this kind of PLAGs are the Galactic bulge 
and rich star fields of the Milky Way spiral arms. However, 
the monitoring of other galaxies, like LMC or M31, should 
give prominence to self-lensing PLAG events. It must be 
said that Protostars, PPNe and other exotic objects are 
not very common in the galaxy, so that their total optical 
depth in observations towards the bulge barely reaches 
10 -11 which is 5 orders of magnitude lower than disk star 
optical depth. 

In a cosmological context, following the idea to find 
extra-galactic MACHOs by monitoring quasars behind the 
Virgo cluster of galaxies (Tadros et al. 1998), or by moni- 
toring M87 with the pixel lensing technique (Gould 1995), 
the search of intracluster PLAG events could give inter- 
esting results. However, considering the actual estimates 
of the abundance of intracluster nebulae (Ciardullo et al. 
1998), the optical depth towards M87 is around 10~ 12 , 
which is still very low. 

5. Summary 

We studied gravitational microlensing due to point-like 
sources surrounded by diffuse gas clouds. Provided the 
size of the gas cloud is about of the order of the corre- 
sponding Einstein radius, we find that the presence of the 
gas cloud, via absorption or refraction effects, can sub- 
stantially modify the classical Paczyhski light curve and 
induce a strong chromaticity. The observation of these ef- 
fects might not be easy, since they require a precise pho- 
tometry and a continuos monitoring, as done for instance 
by the Planet and GMAN collaborations to find planets. 
Otherwise, if one looks for long duration events, the ob- 
servations should last several years, in order to identify 
longer events. 
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We have hypothesized that cold dark PLAGs may be 
present in galactic halos and we have mainly focused on 
the determination of the characteristics they should have 
to be observable. If the existence of PLAGs in the halo is 
proved, then the recent estimates of baryonic dark mat- 
ter should be revised to take into account the presence 
of gas clouds in the halo. Disk PLAG events, caused by 
known nebulae, or events caused by intracluster nebulae 
should be very rare, if the estimates of their abundances 
are correct. 

Due to the chromaticity, the present used selection cri- 
teria for microlensing reject PLAG events. It would thus 
be interesting to check whether already in the existing 
data there are PLAG microlensing events. We believe that 
a search for PLAG type events is necessary to complete 
the picture about baryonic dark matter in the halo. 
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